o 
o 

(N 

a 

(N 



X 



Preprint ESI 1647, TUW-05-07 

Dynamical instability of fluid spheres in the presence of a 

cosmological constant 

C. G. BdhmeiQ 

The Erwin Schrddinger International Institute for Mathematical Physics, 
Boltzmanngasse 9, A-1090 Wien, Austria and 
Institut fiir Theoretische Physik, Technische Universitdt Wien, 
Wiedner Hauptstr. 8-10, A-IO4O Wien, Austria 

T. Harkql 

Department of Physics, The University of Hong Kong, 
Pokfulam Road, Hong Kong SAR, P. R. China 
(Dated: February 7, 2008) 



> 

]0 ■ Abstract 
O 

o 
o 



The equations describing the adiabatic, small radial oscillations of general relativistic stars are 
generalized to include the effects of a cosmological constant. The generalized eigenvalue equation 
for the normal modes is used to study the changes in the stability of the homogeneous sphere 
induced by the presence of the cosmological constant. The variation of the critical adiabatic index 



^P' as a function of the central pressure is studied numerically for different trial functions. The presence 



of a large cosmological constant significantly increases the value of the critical adiabatic index. The 
dynamical stability condition of the homogeneous star in the Schwarzschild-de Sitter geometry is 
obtained and several bounds on the maximum allowable value for a cosmological constant are 
derived from stability considerations. 
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I. INTRODUCTION 



The possibility that the cosmological constant be nonzero and dominates the energy 
density of the universe today is one of the most intriguing problems of the contemporary 
jhysics. Considered already in 1896 by von Seeliger and Neumann (see Schunck and Mielke 
l(] and references therein), who added it to the Poisson equation for the Newtonian potential 
to compensate the energy density of the 'aether', introduced by Einstein in his equations 
of general relativity to obtain a static model of the Universe, eliminated by him after the 
discovery of the redshift of the stars by Hubble, the cosmological constant was reintroduced 
next by Bondi and Gold and Hoyle and Narlikar 4] to resolve an age crisis and to 
construct a universe that satisfied the "Perfect Cosmological Principle". The cosmological 
constant appeared again in inflationary models of the universe introduced by Guth . The 
density of the positive isotropic energy density of the scalar field (infiaton) that dominates 
the first stages of the evolution of the Universe behaves like a cosmological constant leading 
to a rapid cosmological expansion of the universe during a de Sitter phase Gj. 

Within the classical general relativity, the existence of a cosmological constant is equiva- 
lent to the postulate that the total energy momentum tensor of the Universe T^^^^ possesses 
an additional piece T-^\ besides that of its matter content T^^\ of the form T^^^^ = Agik, 
where generally one may assume that the cosmological constant A is a scalar function of 
space and time coordinates. Such a form of the additional piece Tj^^^ has previously been 
obtained in certain field-theoretical models and is interpreted as a vacuum contribution to 



the energy momentum tensor 



Agik = ^ (t^P) = ^ {pv) 9ik, (1) 

where pv is the energy of the vacuum (see Padmanabhan HI] for a recent review of the 
cosmological constant problem). The vacuum value of Tj^ thus appears in the form of a 
cosmological constant A in the gravitational field equations. In the framework of standard 
general relativity the contracted Bianchi identity requires A to be a constant. However, 
generalized physical models with time and space varying cosmological constant (decaying 
vacuum energy) have been intensively investigated in the physical literature [l^ . 

A natural thermal interpretation of the cosmological constant has been proposed based 
on the fact that the de Sitter vacuum can be thought as being hot, since even a geodesic 
observer in this manifold will detect an isotropic background of thermal radiation with 



temperature Ty = (A/127r^)^''^ 13|. The cosmological constant can thus be interpreted as 
a parameter measuring the intrinsic temperature of the empty space-time, or in a sense, of 
the geometry itself. 

The first pressing piece of data, which motivated a reconsideration of the cosmological 
constant involved the study of Type la Supernovae. Observations of Type la Supernovae 
with redshift up to about ^ ~ 1 provided evidence that we may live in a low mass-density 
Universe, with the contribution of the non-relativistic matter (baryonic plus dark) to the 



HQ. 



M Il5|, UJ . The value of 



total energy density of the Universe of the order of f2m ~ 0.3 

is significantly less than unity [itI, and consequently either the Universe is open or there 
is some additional energy density p sufficient to reach the value fitotai = 1, predicted by 
inflationary theory. Observations also show that the deceleration parameter of the Universe 
q is in the range —1 < q < 0, and the present-day Universe undergoes an accelerated 
expansionary evolution ij]. 

Several physical models have been proposed to give a consistent physical interpretation 



to these observational facts. One candidate, and maybe the most convincing one 



or the 



missing energy is the vacuum energy density or the cosmological constant A [2, lUl]. The 
presence of a cosmological constant implies that, for the first time after inflation, in the 
present epoch its role in the dynamics of the Universe becomes dominant jiol . 

On the other hand there is the possibility that scalar fields present in the early universe 
could condense to form the so named boson stars. There are also suggestions that the dark 
matter could be made up of bosonic particles. This bosonic matter would condense through 
some sort of Jeans instability to form compact gravitating objects. A boson star can have a 
mass comparable to that of a neutron star [20]. A detailed analysis of the mass dependence 
of boson stars on the self-interacting potential done by Mielke and Schunck [21!] has shown 
that the mass M of the bosons stars can range from M = 10^^ kg for mini-boson stars 
with radius R ~ 10^^^ m to values that can reach or extend the limiting mass of 3.23M0 
of the neutron stars. For the very light universal axion of the effective string models, the 
gravitational mass is in the range of ~ 0.5Mq. The density of a mini-boson star can exceeds 
by a factor of 10'^^ the density of a neutron star. 

In general, a boson star is a compact, completely regular configuration with structured 
layers due to the anisotropy of scalar matter, an exponentially decreasing 'halo', a critical 
mass inverse proportional to the constituent mass, an effective radius and a large particle 
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number (for a detailed review of the boson star properties see 22|). The simplest kind of 
boson star is made up of a self-interacting complex scalar field $ describing a state of zero 
temperature. The self-consistent coupling of the scalar field to its own gravitational field is 



via the Lagrangian 



22|,|23|,y, 



(2) 



2k " ^" ' 2 

where is the self-interaction potential usually taken in the form V^(|$p) = 

m^|$P + A|$|^/4. Here m is the mass of the scalar field particle (the boson) and A is 
the self-interaction parameter. For the bosonic field the stationarity ansatz is assumed, 
$(t, r) = (j){r) exp(— zcjt). If we suppose that in the star's interior regions and for some field 
configurations $ is a slowly varying function of r, so that it is nearly a constant, then in 
the gravitational field equations the scalar field will play the role of a cosmological constant, 
which could also describe a mixture of ordinary matter and bosonic particles. 

At interplanetary distances, the effect of the cosmological constant could be impercep- 
tible. However, Cardona and Tejeiro [25] have shown that a bound of A can be obtained 
using the values observed from the Mercury's perihelion shift in a Schwarzschild-de Sitter 
space-time, which gives A < 10~^^cm~^. For compact general relativistic objects the cosmo- 
logical constant modifies the mass-radius ratio M/R of compact general relativistic objects. 
Generalized Buchdahl type inequalities for M/ R in the presence of a cosmological constant 



have been derived by Mak et al. 



26| and by Bohmer [27|, see also Stuchh'k 28|]. Thus, for 



a non-zero A, the mass-radius ratio is given by 2M/R < 4/9 + (2/3) (4/9 - AR^/Sy/"^ j27|. 
For a recent discussion of the astrophysical bounds on the cosmological constant obtained 
from the study of spherically symmetric bodies^ee Balagueta-Antohnez et al. Q. 

In a very influential paper, Chandrasekhar [31] studied the infinitesimal adiabatic oscil- 
lations of a gaseous sphere in a general relativistic framework. In particular, from stability 
considerations he derived a condition on the radius R of the stable dense general relativistic 
star, which must obey the constraint R > {2GM /c^) Kj (7 — 4/3), where 7 is the ratio of 
the specific heats of the matter and ii" is a constant which depends on the equation of state 
of the matter. For homogeneous spheres K = 19/42, while for polytropic stars K has values 
in the range K = 0.565 {n = 1) and K = 1.124 {n = 3). The stability against small radial 
oscillations of equilibrium configurations of cold, gravitationally bound states of complex 
scalar fields with equilibrium configurations exhibiting a mass profile against central den- 
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sity similar to that of ordinary neutron stars was studied by Gleiser and Watkins [32|. By 
studying the behavior of the eigenfrequencies of the perturbations for different values of 
the central density a it was shown by using both analytical and numerical methods that 
configurations with central densities greater than a given value crc(O) are unstable against ra- 
dial perturbations. The stability of boson star solutions in a D-dimensional, asymptotically 
anti-de Sitter space-time in the presence of a cosmological term was discussed by Astefanesei 



and Radu 33|]. It was found that for D > 3 the boson star properties are close to those in 
four dimensions with a vanishing cosmological constant. A different behavior was noticed 
for the solutions in the three dimensional case. 

It is the purpose of the present paper to study the effects of the possible existence of 
a cosmological constant on the radial oscillations of a dense general relativistic star. As a 
first step we generalize the pulsation equation to include A. This equation is used to study 
the effect of A on homogeneous compact astrophysical objects. The values of the critical 
adiabatic index 7 are estimated numerically for different values of the cosmological constant 
A. The stability criterion of general relativity for a homogeneous fluid sphere is generalized 
to take into account the presence of the cosmological constant. Based on the analysis of 
the solutions of the gravitational field equations for a constant density sphere some other 
simple criteria of stability for compact objects in the presence of a cosmological constant 
are obtained. 

This paper is organized as follows. In Section 2, using the Einstein gravitational field 
equations, we derive the eigenvalue equation for adiabatic radial pulsations for the case 
A 7^ 0. In Section 3 we discuss the mass continuity and Tolman-Oppenheimer-Volkoff (TOV) 
equations in the presence of a cosmological constant and the solution for the homogeneous 
sphere is obtained. The stability of the homogeneous star is considered in Section 4. The 
results are discussed and summarized in Section 5. 

II. ADIABATIC RADIAL PULSATIONS OF RELATIVISTIC STARS IN THE 
SCHWARZSCHILD-DE SITTER GEOMETRY 

In the present Section we will develop the general formalism for the study of perturba- 
tions of the exact solutions of Einstein's field equations for fluid spheres in the presence of 



a cosmological constant, by following the approach introduced by Chandrasekhar 31|- We 
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shall consider only perturbations that preserve spherical symmetry. Under these perturba- 
tions only radial motions will ensue. The metric of the corresponding space-time can be 
taken in the form 

ds^ = {dx') ' - e^dr^ - [dO^ + sin^ Odcj)^) , (3) 

where v and A are functions of x° = ct and of the radial coordinate r only. 
The energy-momentum tensor for a spherically symmetric space-time is 



(4) 



where e is the energy density of the matter and p is the thermodynamic pressure, = dx^ jds 
is the four-velocity. 

The gravitational field equations corresponding to the line element given by Eq. ([3]) 
and with the cosmological constant included [m], together with the 0th order equations, 
corresponding to the static equilibrium presented in Appendix A . 

Writing A = Aq + 5A, i/ = z/q + bv, p = po + 5p and p = Pq + 6p, where the index refers to 
the unperturbed metric and physical quantities, we find that to the first order the perturbed 
components of the energy-momentum tensor are given by Tq = e, = —p, i = 1,2,3 (no 
summation), Tq = (cq + Po) v and = (cq + Po) ve^°~'^°, where v = dr/dx^. 

By introducing the "Lagrangian displacement" ^ defined as f = d^/dt, we obtain 
the variations of the metric functions and of the energy density in a form similar to 
the A = case: 5\ = - (Af, + u'q) ^, (Su)' = [Sp/ (eo + Po) - (uq + 1/r) ^] (Aq + z/q) and 
6e = — (l/r^) [r^ (eo + Po) ^]', respectively. 

In order to obtain an expression for Sp we need to impose an extra condition on the 
system. This condition is the law of conservation of the baryon number density n, which 
can be written as [nu"^).- = 0. Generally is a function of both energy density and pressure, 
n = n{e,p). By taking n = no{r) + (5n(x°,r) it follows that the variation of the pressure 
can be obtained as 



5p = -^-r 

dr 



IPo 



exp(z/o/2) 



[r2exp(-z/o/2)e]' 



(5) 



where 7 is the adiabatic index defined as [31 1 

dn \ ~^ 
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^ dp) 



n- (e + p) 



dn 
97 



(6) 
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With the use of the previous results for the perturbed quantities, and by assuming that all 
the perturbations have a dependence on x° of the form exp(ia;x°), from the linearized Ein- 
stein field equations we obtain the Sturm- Liouville eigenvalue equation for the eigenmodes 

in the presence of a cosmological constant as 



(7) 



where we have redefined the "Lagrangian displacement" ^ as ^ — > exp[i//2(r)]C (cor- 
responding to a Lagrangian displacement of the radial coordinate of the form 5r(x°, r) = 

r^^ exp[z//2(r)]Cexp(icc;a;°) and we have denoted 



n 



7|g[A(r)+3i/(r)]/2^ 



IP 



P 



/2 



r 



SvrG 



p - A ) {e + p)e 



\{r) 



W = ^ + ^ cf3A(r)+z/(r)]/2_ 



(8) 

(9) 
(10) 



The boundary conditions for C(r) are that C('")/''"^ is finite or zero as r ^ and that the 
Lagrangian variation of the pressure Ap = — (ype'^^^^^'^/r'^) d(/dr vanishes at the surface of 
the star. 

For all stellar models of physical interest the frequency spectrum of the normal radial 
modes is discrete; the n-th normal mode has n nodes between the center and the surface 
of the star. The normal mode eigenfunctions are orthogonal with respect to the weight 
function W{r): WQ(kdr = 0, if j ^ k. 

In the presence of a cosmological constant the system of equations governing infinitesimal 
adiabatic radial pulsations is given by 

C dp 



dr r \ 7p 



e -|- p dr ' 



(11) 



dAp 
dr 



%^e'^ir)-'^^^)^e + p)r-A^ 
& dr 



e -|- p \dr 

1 dp _ AtiG 
e + pdr 



(e -\- p) re 



A(r) 



Ap. 



(12) 
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In Eqs. (fTT!) and (fT2|) ^ is the relative radial displacement ^ = ^r/r^ where Ar is the 
radial displacement. Ap is the corresponding radial displacement of the pressure. Eqs. ffTTj) 
and (IT^ represents the generalization to the case of the Schwarzschild-de Sitter geometry 



of the oscillations equations introduced by Gondek et al. 



35|. 



These pulsation equations have to be integrated with the boundary conditions ( Ap) 

center 

and (Ap)^^rf^^^ = 0. 



center 



III. STATIC FLUID SPHERES WITH NON-ZERO COSMOLOGICAL CON- 
STANT 



The basic equations describing the equilibrium of a static fluid sphere in the presence of 
a cosmological constant are given by Eqs. flA5|) - flA7|) in Appendix A. These equations have 
to be considered together with the equation z/q = — 2]9q/ (cq +Vo)i which can be obtained 
from the condition T-^^ = 0, describing the conservation of the energy momentum tensor. 

With the use of the gravitational field equations and of the conservation equation the 
Tolman-Oppenheimer-Volkoff (TOV) equation, describing the hydrostatic equilibrium of a 



fluid sphere in t 
in the form 26, 



le presence of a cosmological constant A 7^ 0, can be immediately obtained 



29 



36| 



(60+Po) [f M(r) + (^po - t) r^^] 



. 2GM{r) A 2 
J- — 5 — "o"' 



where we have denoted 



M(r) = % f eor^dr. 
c Jo 

For the TOV equation for boson stars see Schunck and Mielke 



(13) 



(14) 



221. 



To obtain a dimensionless form of Eq. (fT3!) and of the mass continuity equation, dM/ dr = 
{An / c^) e^r'^ , we introduce a dimensionless independent variable t] and the dimensionless 
functions eo{ri), Po{ri) and m (r^) by means of the transformations r = ar], eo = eceo^rj), 
Po = ecPoiff) ^"^^ ^ — M*m {■)]), respectively. Here a is a scale factor (a characteristic 
length), ec the central energy density of the star and M* a characteristic mass. With the use 
of these transformations we obtain the mass continuity and TOV equations in the following 
dimensionless form 

dPo (r/) [so [r]) + Po {v)] [^(^) + (^) - 



drj 



7]^ 



2m(r)) 



(15) 
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^ = V'eo iv) , (16) 

where we denoted 

Eqs. (fT3l) and (fT4l) must be integrated with the boundary conditions po (R) = and 
M(0) = 0, respectively, where R is the radius of the fluid sphere, together with an equation 
of state of the form pq = Poi^o)- For a wide class of equations of state one can prove 
uniqueness and existence of these solutions, see 0, Q • 

In the new variables the equation of state becomes Pq = Po{eo), while the boundary 
conditions are given by m(0) = and Pq ('7s') = 0, where rjs = R/a is the value of the 
dimensionless radial coordinate rj at the surface of the star. If we suppose that m is an 
increasing function of 1], while and Pq are decreasing functions of the same argument, 
then it follows that generally Eq G [1, 0] and Pq G [Pc, 0], where Pc = -Po(O) = Po(0)/ec is the 
value of the pressure at the center of the star. 

A particularly important fluid sphere configuration is the one corresponding to the con- 
stant density case, e = ec =constant. This condition gives Eq = 1, Vr/ G [0,r7s'], and the mass 
continuity equation can be integrated immediately to give 

3 

m{v) = y. (18) 

The solution of the TOV equation with the boundary condition Pq (0) = Pc is given by 

_ (l/3^0a-Vl-(2/3 + 0„^ 

where we denoted 

a = , ^ . (20) 

1/3-1 + Pc ^ ^ 

The radius R of the star can be obtained from the condition Pq (ris) = and is given by 

„^ V6Pc{2Pc + l)-9Pc{Pc + 2)l 

^ 3P, + 1-3/ ■ ^ ^ 

The values of the central pressure Pc depend on the physically allowed upper limit for the 
pressure. If we consider the classical restriction of general relativity, p < e/3, asp > 0, it 
follows that Pc G [0, 1/3]. A more general restriction can be obtained by assuming that p < e, 
the upper limit being the strong equation of state for the hot nucleonic gas. It is believed 



9 



that matter actually behaves in this manner at densities above about ten times the nuclear 
density, that is, at densities greater than 10^^ g/cm^ and at temperatures T = (e/a)^^^ > 10^^ 
K, where a is the radiation constant In this case Pc € [0, 1]. Hence generally < 
where w e [0, 1]. Consequently, a e [3/(1 - 3/), 6/(4 - 3/)]. 
The metric coefficient exp (— Aq) is given by 



exp [-Ao (r^)] = 1 - (2/3 + /) r,'. 



(22) 



From the Bianchi identity and the matching across the vacuum boundary of the star it 



follows that exp [u [r])] = [1 — (2/3 + I) r/l] / [1 + Pq (r/)]^, or, equivalently. 



(2/3 + ly 



e 2 



1 2 



a 



(23) 



IV. STABILITY OF HOMOGENEOUS STARS IN THE PRESENCE OF A COS- 
MOLOGICAL CONSTANT 



The eigenvalue problem formulated in Eq. ([7]) can be re-expressed in a well-known vari- 
ational form: the extremal values of the quantity 



i: n(f)^-gc^ 



dr 



(24) 



are the eigenvalues uj"^ of Eq. (JTj) and the functions C,{r) which give the extremal values are 
the corresponding eigenfunctions. A sufficient condition for the dynamical instability of a 



mass is that the right hand side o: 
satisfies the boundary conditions 



Eq. ([2^ vanishes for some chosen trial function C, which 

3il, m 



In terms of the system of variables introduced in the previous Section Eq. f l24l) takes the 
form 

"VS 



vs 



2 2 

a uj 



e 2 



1 + Po a^^2, 

V 



\+3v 

V 



7]^ 



p/2 

-'^0 



i^_2('Po--/)(l + Po)e^ 



(25) 



Setting a;^ = and assuming that the adiabatic index is constant throughout the star, 
we obtain the value of the critical adiabatic index 7c as 



/c 



VS e^Z~ 




7c 



p/2 

1+Po 



^-2(Po-|0(l + ^o)e^ Cdv 



(26) 
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For 7 < 7c dynamical instability occurs and the star collapses. For a homogeneous star in 
the presence of a cosmological constant the equilibrium metric coefficients and the pressure 
distribution have been obtained in the previous Section. In order to find the range of values 
of 7c we need to chose some explicit functional forms for (. For the trial function ( = the 
variation of 7c as a function of the central pressure Pc is represented, for different values of 
I, in Fig. 1. 



25 
20 
15 
10 
5 




0.2 0.4 0.6 0. 



FIG. 1: Critical adiabatic index 7c as a function of the central pressure Pc for the trial function 
( = r]^ and for different values of the cosmological constant: I = (solid curve), I = 0.1 (dotted 
curve), / = 0.2 (dashed curve) and I = 0.3 (long dashed curve). 



In terms of the gravitational radius of the star Rg = 2GMs/c^, where Ms = 4:7recR'^/c^ 
is the total mass of the star, and of the radius R, the critical adiabatic index 7c can be 
represented as a power series in Rg/ R of the form 



7c 



^-l 19 
1 - 3/ ^ 42 



21, 

1 / 

19 



Rg 

R 



146 
^441 

+0 



. 87,/ 21, 

1 n iH — / 

73 V 58 



R„ 



Rg_ 

R 



R 



(27) 



For the trial function C, — rj^ exp{—u / A) Fig. 2 represents 7c. 
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FIG. 2: Critical adiabatic index 7c as a function of the central pressure Pc for the trial function 
C, = r/^ exp(— and for different values of the cosmological constant: / = (solid curve), / = 0.1 
(dotted curve), / = 0.2 (dashed curve) and / = 0.3 (long dashed curve). 

In this case the representation of 7c in terms of Rg/R is given by 



19 



7c 



+ 



3/ ^ 42 



5401 



21, 

1 / 

19 



15876 



26691 
21604" 



I I 



Rg 

R 



447 
1271' 



Rg_ 

R 



+0 



Rg 

R 



(28) 



In first order in Rg/R both Ens. ( l27j) and ( l28i) give the same result. This shows that the 
Chandrasekhar stability limit 3l| is generally independent of the form of the trial function. 

Lastly, 7c as a function of both / and Pc and for the trial function ( = rj^, is represented 
in Fig. 3. 



For a homogeneous sphere and for the trial function ( = f]^ all the integrals in Eq. ( !25|l 
can be calculated exactly. The corresponding pulsation equation and the details of the 
calculations are presented in Appendix B. 

Performing a series expansion of Eq. fIBll) in Appendix B, keeping only the terms up to 
the second order in ^5 and neglecting the terms containing P we obtain the equation 



4 390^ _ 406 -3 (142 + 3037)/ ^2 

7 37 — 1 / = 775 T Oq. 

' 3 ^ ' ' 378 (1 + /) ^ 



Taking / = gives 



7 



4 _ 39O7-4O6 
3 ~ 378" 



-Vs- 



(29) 



(30) 
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FIG. 3: Critical adiabatic index 7c as a function of the central pressure Pc and of I for the trial 
function C, = rf 



Therefore 

4 19 2GMs 
3 ^ 42 ^ c^R 

or, equivalent ly, 



7-o<^x--/, (31) 



19 

R < — . (32) 

42(7-4/3) ^ ^ ' 

Hence from Eq. fl29|l we have recovered the stability condition derived by Chandrasekhar 



3l| for A = 0. For Z 7^ we obtain 



o ^ 19 1 - 273//19 

^ < 427-4(l + 9//4)/3^^- ^^^^ 

To obtain Eq. (133|) we have approximated 7 in all terms containing the product 7/ by 
7 = 4/3. Eq. ( l33l) represents the condition for the dynamical stability of a massive general 
relativistic object in the presence of a cosmological constant. It generalizes to the case A 7^ 



the relation initially derived by Chandrasekhar 31[. In the limit / = we obtain again Eq. 

Alternatively, we may divide Eq. (!25|) by the integral on its left-hand side and proceed 
as described above. Up to terms of the order ^| and / this yields 

4 547-53- 15(4 + 97)/ ,2 
=^-3-^(37-1) 63(2/3 + /) (34) 

which generalizes the original Eq. (77) of Chandrasekhar [31]- 
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Approximating, as before, 7 in all terms containing 7/ by 7 = 4/3 and letting u we 
arrive at 

19 1 - 240//19 
^ < 427-4(l + 9//4)/3^^- ^^^^ 

For the trial function ^ = ?7'^exp (— z//4) the integrals in Eq. (123]) can also be calculated 
exactly. A series expansion in terms of 6s of the result for this choice gives 

7 - i - (37 - 1) i = 1237- 107- 3(65 + 1027)/^,^ 

' 3 ^ ' ^ 189(2/3 + /) ^ ^ 

For / = Eq. (136|) gives the same limit as obtained for the case of the trial function 
C, = rj^, represented by Eq. fl32|) . For / 7^ 0, by approximating again 7 by 4/3 in the terms 
containing the product 7/, we obtain 

^ ^ 427-4(l + 9//4) /3^^' ^^^^ 
which differs only slightly in the numerical value of the coefficient of / from the previous 
result. As was already noted, the slight differences in the approximations show that firstly, 
the stability limits are generally independent of the trial function and that secondly, they 
are quite unaffected by the 'location' where the higher order terms are neglected. 

By imposing the condition 7 < 7^ in Eqs. (1271) and fl28|) gives an other form of the stability 
criterion for high density relativistic stars in the presence of the cosmological constant 

„ 19 1-21//19 

R < tR„. 38) 

42 7 -(4/3 -0(1 -3/)"' ' ^ ' 

The denominator of the latter stability criterion coincides with those derived before, if 
terms containing are neglected in a series expansion. 



V. DISCUSSIONS AND FINAL REMARKS 

The hypothesis that a non-zero cosmological constant could play an important role in 
the structure of compact stellar objects seems to be difficult to accommodate with the 
smallness of the present determined or postulated values of this physical quantity. While 
a non- vanishing cosmological constant, if definitely confirmed, would carry a significant 
implication for our understanding of the structure and global dynamics of the Universe as a 
whole, the ways in which one can clearly test it in an astrophysical or astronomical setting 



are very limited. Peebles |3y] has shown that the cosmological constant has little effect 
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on the local dynamics. The effects on the galactic dynamics of A can hardly be observed 
also with distant clusters, because it appears to be always cancelled in observable quantities 



40(1 . But, on the other hand, the effects of the cosmological constant could be important 



for stellar objects in the far past, most of the cosmological models with decaying vacuum 
energy predicting a time variable and decreasing cosmological constant 411]. Therefore, the 
interior structure of early stars could have been strongly influenced by the presence of a 
cosmological constant. 

In the presence of a cosmological constant the radius R of the homogeneous stellar con- 
figuration is related to the gravitational (Schwarzschild) radius of the star Rg by the relation 

o_ (3Pe + l-3/f 

P,[2(2P, + l)-3(Pe + 2)/]'^^- ^^^^ 

From the restriction Pc < 1 on the pressure we find the following simple stability criterion 
for a homogeneous star in the presence of a cosmological constant 

R > ,^,Rr (40) 

- w[2{2w + l)-3{w + 2)l] ^ ^ ^ 

For w = 1/3 we have P > 9 (2 - 3/)^ Rg/ (10 - 21/), while for w = 1 we obtain R > 
{A-3lfRg/[3 (2-3/)]. 

The condition of the non-negativity of the radius of the star, P > 0, imposes the constraint 
/ = Ac^/12'kG€c < 2{2w + l)/3{w + 2) on the cosmological constant, or, equivalently, 

A < (!!^?^. (41) 

For a star composed from matter having a density equivalent to the typical nuclear 
density, 2 x lO^^g/cm^ Eq. dUl) leads to an upper limit of the cosmological constant equal 
to 

A < 3 X lO^^^cm-^ (42) 

In the presence of a cosmological constant the dynamical stability condition is modified 
according to Eq. (!33l) . The critical value 7^ = 4/3 of the adiabatic index is replaced by 
7^ = (4/3)(l + 9//4) = (1 + 3Ac'^/167rG'ec). A large value of / can increase significantly the 
value of 7^. This could be the case for mixtures of the potential energy dominated bosonic 
matter and usual matter, when the cosmological constant A can be interpreted as the static 
equilibrium value of the scalar field potential, A = \^(]$ol^) = const. In this case / could have 
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a large value, even of the order of unity. The presence of the bosonic component strongly 
affects the dynamic stability of the system. The critical values of 7c of configurations with 
given ec require that the radius of the star exceeds a certain lower limit i?, if the configuration 
is to be stable. In the presence of a cosmological constant instability will occur when 7 is in 
excess of 4(1 + 9//4)/3 by a small amount. Since for a perfect gas the maximum admissible 
value for 7 is 7 = 5/3, the stability condition imposes the restriction 

9^ 



5-4(1+4 



> 0, (43) 



or, equivalent ly, 

AnGec , , 

A < (44) 

By using the same restriction on 7 in Eq. fl38l) we obtain 

A < (45) 

These conditions, obtained from stability considerations, differs from the condition given 
by Eq. fj4T|) only by factors of the order of unity. 

The properties of boson stars and the influence of the cosmological term on their proper- 
ties in arbitrary -D-dimensional, asymptotically anti de Sitter space-times have been recently 



considered in 33|] . By assuming for the spherically symmetric bound state of the scalar field a 
stationary ansatz $ = 0(r) exp {—iat) the metric for the boson star can be chosen in the form 
ds'^ = F{r)dr'^ + r'^dVt\,_2 — F{r) exp [—25{r)] dt^, where dVt\)_2 = 0Jabdx"-dx^ is the line ele- 
ment of a unit {D — 2)-dimensional sphere and F{r) = l—2'm{r)/r^~^—2Ar'^/{D—2){D~l). 
The function m(r) is related to the local mass-energy density up to some dimension- 
dependent factor. With the help of these parameters one can define the particle number 
N = {47r(^-i)/Vr [{D - 1) /2]} (j)^e^F-^r^-^dr and the effective radius of the boson 
star as R = {47r(^-i)/Vr [{D - 1) /2] A^} 4>'^e^F-^r^~^dr. For A = 0, boson stars are 
characterized by an exponential decay of the scalar field, for which the mass term in the 
potential is responsible. For A 7^ the situation is quite different; for r 00, one finds that 
m{r) ~ M + Ar'^'-''~^^~^ , where A and C are some constants depending on the dimensionality 



of the space-time D and A, whose exact form is given in [33|], and M is the ADM mass of the 
star. Thus, the cosmological constant implies a complicated power decay at infinity rather 
than an exponential one as found in the asymptotically flat space case. The maximum mass 
^max and the maximum particle number A^max for the boson star decreases with the value 
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of the cosmological constant, but the general properties of the solution are the same as in 
the A = case. 

The analysis of the stability of the boson stars follows the standard approach, by assuming 
that the scalar field and the metric tensor components can be perturbed in the usual way, 
with the perturbation of the scalar field written as 50 = fi{t, r)+i(j)Q{r)g (t, r). Then the field 
equations and the Klein-Gordon equation can be linearized to a system of two self-adjoint 
coupled equations, the pulsation equations for the boson stars. These equations have to be 
integrated with the boundary conditions r'^g' — at the origin and /i ^ 0, g'r'^'^+^ 
for r — s> oo. 

By assuming that all perturbations are of the form exp {ixt) , where x is the characteristic 
frequency to be determined, one can again formulate the stability problem of the boson stars 
in the presence of a cosmological constant as a variational problem, with the eigenvalues 
forming a discrete sequence Xo — Xi ^ ••• ^ Xn- The numerical results show that there is a 
critical value of the central density of the star 0c(O) so that > for central values of the 
scalar field smaller than 0c(O) while other configurations are unstable. 

Due to the close mathematical and even physical analogy between the two eigenvalue 
problems for boson and "normal" stars, we expect that stability conditions of the form given 
by Eqs. ( 133!) . (135!) and (!38|) should also exist in the case of boson stars. For boson stars the 
corresponding mass appearing in the stability conditions is the ADM mass M, and the radius 
is defined as above. The coefficient 7 becomes a scalar field-dependent quantity 7^5, which 
could be given as a function of the the central density of the star 0c(O), 7_bs = 7 [0c(O)]. The 
values of the numerical coefficients could be also very different as compared to normal stars. 
Therefore for boson stars a mass-radius relation of the form R < [f{l)/{'^BS ~ const)]M 
is a direct consequence of the variational formulation of the stability problem and of the 
boundary conditions, requiring the vanishing of some quantities at center of the star and 
infinity, respectively. But finding exact analytical stability conditions for boson stars requires 
further mathematical and numerical investigations. 
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APPENDIX A: FIELD EQUATIONS 

For a spherically symmetric matter distribution the gravitational field equations are given 

^re-^y ^l_^Ty-Ar\ (Al) 

, e^-1 SnG. , , 
u' = — —T^r - At, (A2) 

1 .f„ u'^ u' - X' u'X' 



2 V 2 r 2 



-A 



A SttG 



e — = 



Tl (A4) 



The zeroth order (or static equilibrium) equations are 

(A5) 

Cr 

ly'o = + i-^Po - Aj r, (A6) 

where the index refers to the unperturbed metric and physical quantities. 



APPENDIX B: PULSATION EQUATION 

In this Appendix we present the pulsation equation for a homogeneous sphere and for 
the trial function ( = rf. By introducing the new variable 6 = arcsin ( a/2/3 + /r^) , denoting 



6s — arcsin (a/2/3 + Irjs) and taking into account that a = cos6'5/(l/3 — /), in the limit 
a; — > we obtain the following equation for 9s'- 

360 [22 - 347 + 9/ (1 + 147 - 9/7 + 3 (2 + 7) f)] 9s cos (^5) + 
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1080 [2 + 31 + 37 (3Z - 1)] 9s cos (8^5) - 
20 [134 - 63/ + 594/2 27 (3/ - 1) (7 - 8/ + 6/^) 7] sin {9s) + 
5 [-440 + 7177 + 9/ (-22 (2 + 3/^) + 9 (-31 + / (13 + /) 7))] x 

sin (3^5)- 

[-10 (3/ - 4) (2 + 3/)^ + 27 (3/ - 1) (7 + (8 + 3/) /) 7] x 

sin (5^5) = 0. (Bl) 
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